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In this paper it is shown that an i(p'^ field theory is a physically acceptable field theory model 
(the spectrum is positive and the theory is unitary). The demonstration rests on the perturbative 
construction of a linear operator C, which is needed to define the Hilbert space inner product. The 
C operator is a new, time-independent observable in PT-symmetric quantum field theory. 
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A scalar gip^ field theory is often used as a pedagogi- 
cal example of perturbative renormalization even though 
this model is not physically realistic (the energy is not 
bounded below). However, we argue that when g = ie 
is imaginary, one obtains a fully acceptable quantum 
field theory and we show how to construct perturba- 
tively the Hilbert space in which cubic scalar field the- 
ories in (D + l)-dimensional Minkowski space-time are 
self-adjoint. Consequently, such theories have positive 
spectra and exhibit unitary time evolution. We consider 
here such complex field-theoretic Hamiltonians as 



H 



12 2 
2f^ '/'x 



(1) 



where J ^ J d^x. [We suppress the time variable 
t in the fields and write Lp{'x.,t) = (px-] The fields in 
(HJ satisfy the ETCR [(^(x, i), 7r(y, t)] = i5{^-y). As 
in quantum mechanics, where the operators x and p 
change sign under parity refiection V, we assume that 
the fields are pseudoscalars and also change sign under 
V: Vip{yi,t)V = -(p(-x,i) and -P7r(x,<)-p = -7r(-x,i). 
Field-theory models like that in Q are of physical inter- 
est because they arise in Reggeon field theory and in the 
study of the Lee- Yang edge singularity . 

This paper is motivated by the observation that the 
cubic complex quantum-mechanical Hamiltonian 



i^^x^ -I- iex^ (e real) 



(2) 

has a positive real spectrum |^ |M Although this 
Hamiltonian is not Hcrmitian in the conventional sense, 
where Hermitian adjoint means complex conjugate and 
transpose, it still defines a unitary theory of quantum 
mechanics This is because in (O is self-adjoint 

with respect to a new inner product that is distinct from 
the inner product of ordinary quantum mechanics '^l. 

For quantum-mechanical theories like that in two 
issues need to be addressed. First, the spectrum of H 
must be shown to be real and positive. In Ref. it 
was observed that spectral positivity of H is associated 
with unbroken space-time reflection symmetry [VT sym- 
metry). [The term unbroken VT symmetry means that 
every eigenstate of H is also an eigenstate of VT . This 
condition guarantees that the eigenvalues of H are real; 

in (O has an unbroken VT symmetry for all real e.] 



Second, it is necessary to construct a Hilbert space with a 
positive norm on which the Hamiltonian acts and to show 
that H is self-adjoint with respect to this nonstandard 
inner product. Using a coordinate-space representation, 
this construction was made in Ref. [5j by introducing a 
new linear operator C: 



C{x,y) 



(l)„{x)(t>niy), 



(3) 



where (j)n{x) {n = 0, 1, 2, . . .) are the eigenstates of H. 

For the cubic Hamiltonian in Q, (t>n{x) satisfies the 
Schrodinger equation 

-i<(a;) -I- i/i2a;20„(a;) -I- iex^(t)n{x) = £;„0„(x) (4) 

and the boundary conditions (j)n{x) — > as a; — > ±oo. 
The eigenstates of H are also eigenstates of VT and are 
normalized so that VT(j)n{x) — (/)*(— x) = (j)n{x). The 
inner product involves CVT conjugation 



CVT 



dxij^'^^{x)x{x), 



where V'''^^(a;) = J dyC{x,y)^p*{—y), instead of conven- 
tional Hermitian conjugation, where the inner product is 
(V'lx) = J dx {x)xix) ■ The novelty of these complex 
Hamiltonians is that the Hilbert space inner product is 
not prespecified; rather, it is dynamically determined by 
H . This new kind of quantum theory is a sort of "boot- 
strap" theory because one must solve for the eigenstates 
of H before knowing what the Hilbert space and the as- 
sociated inner product of the theory are. 

The key problem in understanding a complex Hamilto- 
nian like that in or ||2Jl is to determine the operator C. 
In Ref. perturbative methods were used to calculate 
C to third order in e for the H in The procedure was 
first to solve the Schrodinger equation Q for (/>„ (x) as a 
series in powers of e. Then, 0,i(x) was substituted into 
and the summation over n was performed to obtain 
C{x,y) to order e'^. The result was complicated, but it 
simplified when C was rewritten in exponential form: 



Cix,y) = (e 



S{x + y)+Oie''), 



(5) 



where Q2n+i{x,p) are differential operators depending on 



X and p - 



dx 



In the representation (jSJ only odd pow- 



ers of e appear in the exponent, the coefficients are real. 
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and the derivative operators act on the parity operator 
V ^5{x + y). Also, CV = 3"'^^+''^'^^+"' is Hermitian. 

Calculating C by direct evaluation of the sum in is 
difficult in quantum mechanics because it is necessary to 
determine all the eigenfunctions of H. In quantum field 
theory such a procedure is impossible because there is no 
analog of the Schrodinger eigenvalue problem Q). 

The breakthrough reported here is the discovery of a 
new and powerful method for calculating C by seeking an 
operator representation of it in the form C = e^^^^'^^V. 
This representation, which is suggested by ((SJ), has the 
advantage that Q{x,p) is determined by elementary op- 
erator equations so that the eigenfunctions of H are not 
needed to find Q. Thus, the technique introduced here 
generalizes to quantum field theory. 

We illustrate the generality of the representation C = 
e^V by using two elementary examples: 

Example 1: The complex 2x2 Hamiltonian 



H 



re 
s 



(6) 



taken from Ref. 0, is VT symmetric, where 'P — (il) 
and T is complex conjugation. A nice way to express C 
in the region > sin^ 9 of unbroken VT symmetry is 

C = e'^V, Q = 5f72ln[(l -sina)/(l +sinQ;)] , 



where sin a — (r/s) sin 9 and <72 — q] 
matrix. As 9 0, H becomes Hermitian and C 



Example 2: The Hamiltonian H = ^p^ - 



2^ 



is the Pauli 
- iex has 

an unbroken VT symmetry for all real e. Its eigenvalues 
En — n + ^ + are all real. The C operator is given 
exactly by C = e'^V, where Q = —ep. Again, in the limit 
e ^ the Hamiltonian becomes Hermitian and C ^ V . 

Our new procedure for constructing C for a given H 
is based on three observations made in Ref. [5j: (i) C 
commutes with the space-time reflection operator VT ^ 
but not with V or T separately; (ii) the square of C is 
the identity; (iii) C commutes with H. To summarize. 



For all cubic Hamiltonians, Q{x,p) may be expanded 
as a series in odd powers of e as in ((SJ: Q{x,p) = 
eQi{x,p)A-e'Q3{x,p)+e^Q^{x,p)+- ■ ■ . (In quantum field 
theory we can interpret the coefficients Q2n+i as inter- 
action vertices of 2n -I- 3 factors of the quantum fields.) 
Substituting this expansion into © and collecting the 
coefficients of like powers of e, we obtain a sequence of 
operator equations that can be solved systematically for 
the operator-valued functions Qn{x,p) {n = 1,3,5,...) 
subject to the constraints that Q{x,p) be even in x and 
odd in p. It is interesting that even powers of e yield re- 
dundant information because the equations arising from 
the coefficient of e^" can be derived from the equations 
arising from the coefficients of e^"~^, e^"~'^, . . . , e. The 
first three operator equations for Qn are 



-2ffi, 



-^([Qi,[Q3,i?i]] + [g3,[Qi,i?i]] 



(9) 



For in ||2Jl we solve the equations in lO by substi- 
tuting the most general polynomial form for Qn using 
arbitrary coefficients and then solving for these coeffi- 
cients. For example, to solve the first equation in 0, 
[Hq, Qi] — —2ix^, we take as an ansatz for Qi the most 
general Hermitian cubic polynomial that is even in x and 
odd in p: Qi{x,p) = Mp^ + Nxpx, where M and N are 
unknown coefhcients. The operator equation for Qi is 
satisfied if M = and N = -2^-2. 

To present the solutions for Qn{x,p) we use Sm,n to 
represent the totally symmetrized sum over terms con- 
taining m factors of p and n factors of a; 0. Thus, 



So,o 



1, So,; 



Si, 



Uxp + px). Si:, 



■| {x^p + xpx + px^) , and so on. We have solved Q for 
Qij Qsi Qbj and in closed form in terms of Sm,n'- 



Qi 



128,-10^5 I 40, -8 ( 



(i) [C, VT] = 0, (n) = 1, (in) [C, H] = 0. (7) ^3 = ^A^^^V + f A^-«^3,2 + 8^i-^S^., 12^,-^p, 



Substituting C — e^V into condition (i) in (0), we ob- 
tain e'3(="'P) = 'pTe^^'^ P^VT = e'^^-^-P), so Q{x,p) is 
an even function of x. Next, we substitute C = e^V 
into condition (u) in and get e'5(^.P)-Pe'5(^.P)'P = 
^Q{'j:,p)^q{-x,^p) ^ implies that Qix.p) = 

—Q{—x,-'p). Since Q{x,p) is an even function of x, it 
must also be an odd function of p. Finally, substitut- 
ing C = e'3(^'P)'P into condition (iii) in O, we obtain 
eQ(^,rt[7?,i/] + [e'3(^'P),i/]75 = 0. AU of the Hamiltonians 
El considered in this paper have the form iJ — Hq + cHi , 
where Hq is a free Hamiltonian that commutes with the 
parity operator V and Hi represents a cubic interaction 
term, which anticommutes with V. Hence, for quantum 
theories described by cubic Hamiltonians, we have 



Q 



320,-16 7 544, -14c 512,,-12q 
fi/l,,-10<? 24 736 , .-14^3 , 6368 ,, -12 c 



3,4 



l,,-18c 
-/i b; 



553984 ,,-22„9 , 9 7 792 ,.-20 c , 377344, ,-18, 
35 

315 ^ ""'^3,6 + nf^M ^'*'5'l,8 



Qt - -3T5-M P + 35 

721 024 ,,-16 g_ „ _|_ 1 792 ,,-14 



^ ' 105 
2209024 ,,-20 „5 
105 ^ P 



5,4 



2875648 ,, -18 c 390 336 ,,-16 c , 46 976 ,,-18,„ n r\\ 
105 M -^3,2 35 — M Jl,4 H 5 — M P- (,tUj 



(8) 



These results constitute a seventh-order perturbative ex- 
pansion of C in terms of the operators x and p. 

The representation C ~ e'^V bears a strong resem- 
blance to the WKB ansatz, which is also an exponential 
of a power series. Like our calculation of C here, WKB 
methods determine the energy eigenvalues En to all or- 
ders in powers of h via a system of equations like those 
in ^ without ever using the eigenfunctions (t)n{x) |lOl |. 
Moreover, only the even terms in the WKB series are 



3 



needed to find i?„; the odd terms in the series are re- 
dundant and provide no information about En ^3 • The 
difference between a conventional WKB series and the 
series representation for Q(x,p) is that the first term in 
a WKB series is proportional to h^^ while the expansion 
for Q contains only positive powers of e. However, based 
on the nonperturbative calculations in Ref. Q , we believe 
that for a T^T-symmetric —ex^ theory, the first term in 
the expansion of Q will also be proportional to e^^. 
The massless (strong-coupling) limit /i of C for 
in (^3) is interesting because it is singular. [The di- 



so nega- 



mensionless perturbation parameter is e/i"'"^/^, 
five powers of /i appear in every order in pOfl . Thus, 
as fi 0, the perturbation series for Q ceases to exist.] 
As — > 0, a new nonpolynomial representation for Q 
emerges. To find C when /i = we return to the op- 
erator equations in J^J and seek new solutions for the 
special case -ffo ~ \p^- The situation here is like that 
in Ref. 0, where the objective was to calculate the time 
operator in quantum mechanics. Here, as in the case 
of the time operator, we use Weyl ordering to general- 
ize Sm,n from positive m to negative m. For example, 
S-i^i — \ {x\ + ix). The solution to the first equation 



in ©, [\p^,Qi\ = -2ia;3, is Qi \S-i^i + aS-^fi, 
where a is arbitrary. This solution is odd in p, even in 
x, and has the same dimensions as Qi in (|10|l . Also, 

n — f 3 _L 305^^ c (-135 I 5773^ 75^2^ o , 

^3 - [32+— Oi) ^-llA - lT6- + -8-"-T2" ) -^-13,2 + 

(^+20a) 5-9,6 - -§2^-7^s + -^3-5^10 + pS^i5fl, where 
/3 is a new arbitrary number. 

The new operator techniques introduced in this paper 
extend to systems having several dynamical degrees of 
freedom. We have calculated C to third order in e for 

H = ^{pl+pl) + '^{x^+y^)+iex^y, (11) 

which has two degrees of freedom. The result is 



405 



405 



1 ^1 r,\ — 512^,2^3 I 512 4„ I 1088 c rp 

Jlrp , 512 
■P ^1.2 + 405 

i^2rp 256 q „a I e08 c „ 128^4 8„ /-in-. 



256 „2^ I 512 Q „ I 288 Q „ 32 „2 S 
-405P Jl,2 + 405*^3,1^ + 405*^2,29 - 405a: q 



736 „2r 



where Tm,n represents a totally symmetric product of m 
factors of q and n factors of y. For the Hamiltonian 

H=^{pI+pI+pI) + ^ {x^ +y^ + z') + iexyz, (13) 

which has three degrees of freedom, we have 

Qi{x,y,z,p,q,r) = -^{yzp xzq + xyr) - ^pqr, 
Q3{x,y,z,p,q,r) = ^ {p^ qr + q^ pr qp) 
-^^[pxp{yr + zq) + qyq{xr -|- zp) -|- rzr{xq + yp)] 
--^{xpxqr yqypr + zrzpq) -|- ^{xpxyz -f yqyxz 
-\-zrzxy) — ^ [x^{yr + zq) + y^{xr + zp) 
+z^{xq + yp)] - {p^yz + q^xz + r^xy) . (14) 



Note that the technique used in Ref. to calculate 
C gives Qi in ((T^ and (fT^ . but it becomes hopelessly 
difhcult beyond first order because one encounters the 
problems associated with degenerate energy levels. Our 
new method works to all orders of perturbation theory 
because the eigenf unctions are not needed to obtain C. 

We now show how to apply the powerful techniques we 
have developed in quantum mechanics to the calculation 
of C in quantum field theory. As before, we express C 
in the form C = e^'^i+e Qs+.-.-p^ where now Q2n+i {n ~ 
0, 1, 2, . . .) are real functionals of the fields (^x and tTx. 
To find Qn for iJ in Q we must again solve the operator 
equations in |(5J). In terms of the inverse Green's function 
= (^^ - Vx)(5(x - y) the first equation in Q is 



= -Aijdxipl. (15) 



dx7r,, + yycfxdy<y9xGxy<y3y,(5i 
Recalling the result for Qi in (|10|) we try the ansatz 
1 = dxdydz {Mi^^y^)T:^iTyT:^ + iVx(yz)</'y'^x<;5z) ■ 



The notation -M(xyz) indicates that this function is to- 
tally symmetric in its three arguments and the notation 
A'x(yz) indicates that this function is symmetric under 
the interchange of the second and third arguments. The 
unknown functions M and A'' are form factors; they de- 
scribe the spatial distribution of three-point interactions 
of the field variables in Qi. The nonlocal spatial interac- 
tion of the fields is an intrinsic property of C. 

To determine M and N we substitute Qi into (|15(l and 
obtain a coupled system of differential equations: 

(/i^ - V2)A^x(yz) + (/^' - Vy2)Ay(xz) + (m' - V2)A^,(xy) 

= -6(5(x-y)(5(x-z), 
3(m' - V2)M, 



x(yz) 



A, 



y(xz) 



'(wxy) 



(16) 



We solve these differential equations by Fourier trans- 
forming to momentum space, to obtain 



M, 



(xyz) 



(277) 



2D 



dp dq- 



j(x-y)-p+i(x-z)-q 



^(P,q) 



where D(p,q) = 4[p^q^ — (p • q)^] + 4/x^(p^ 
q^) + 3^^ is positive, and 



A". 



(y^) 



3 (Vy • V, + \^?) Af(xyz) 



(17) 

p q + 

(18) 



For the special case of a (1 + l)-dimensional quantum 
field theory the integral in 117|l evaluates to 



Mr. 



3^^* ) Ko(Aii?), 



(19) 

?2 



'(xyz) — 

where Ko is the associated Bessel function and Fi^ = 
i[(x-y)2 + (y-z)2 + (z-x)2]. 

We can perform the same calculations for cubic quan- 
tum field theories having several interacting scalar fields. 
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Consider first the case of two scalar fields ip^'' and i^x^-* 
whose interaction is governed hy H = H^^'' + h'^^ + 

ie J c?x [</Jx^'']^</'x^\ which is the analog of the quantum- 
mechanical theory described by H in Here, 



To determine C to order e we introduce the ansatz 



(2) 



,A^(2) 7.(^)09(1)^(1) + A/ , (2) (1) (1) 



where Mx(yz), A^^y'z, and ivi^^ 



x(yz) 



are unknown functions 



and the parentheses indicate symmetrization. We get 
4 



x(yz) 



{2n 



\2D 



gi(x-y)-q+i(x-z)-r 



where 2?(q, r) = 4[q2r2 - (q • r)^] + 4/i^(q + r)^ - ■ 
r — fi2+ 4^i^2' ^'^d 



AT 



xy 

(2) 

x(yz) 



Ml 



1 ,,2 



Ai2j ^4(yz)- 



(21) 



For three interacting scalar fields whose dynamics is de- 
scribed hyH = H^g'^'>+ 7?^^) + i?^^) +iejdx ip^^'' tp'i'^ tp^^^ , 
which is the analog of H in we make the ansatz 

+iVg),4')^W^(2) + Mxyzvr(i)4^)^(3) 

The solutions for the unknown functions are as follows: 
-^xyz is given by the integral (|17|l with the more general 
formula /^(p, q) = 4[p2q2 - (p • q)^] -f 4[^^(q2 -|- p • q) -|- 
mKp^ + P ■ q) - M3P ■ q] + with = 2iilnl + 2^2^§ 4- 
2/^2/^3 ~ lA~ lA~ Ms- The coefficients are expressed 
as derivatives acting on M: 



xyz 
xyz 
xyz 



[4Vy • V, 



2(M2 



\+t4- 

[-4Vy • Vz - 4V2 + 2(^2 ^ ^2 _ ^2)j 



'4Vy • - 4V2 + 2(/x2 + ^2 _ ^2)j 



Our new perturbative calculation of C for cubic quan- 
tum field theories is an important step in our ongoing 
program to obtain new physical models by extending 
quantum mechanics and quantum field theory into the 
complex domain. The operator C is a new conserved 
quantity in quantum field theory and this operator is re- 
quired to construct observables and to evaluate matrix 
elements of field operators. 

We hope to generalize the breakthrough reported in 
this paper to noncubic T'T-symmetric quantum field the- 
ories such as a (3-|-l)-diniensional —g^p^ field theory. This 



remarkable model field theory has a positive spectrum, 
is renormalizable, is asymptotically free |l2l |. and has a 
nonzero one-point Green's function Gi = {(f). Conse- 
quently, this model may ultimately help to elucidate the 
dynamics of the Higgs sector of the standard model. 
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